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EXECUTIVE SUMMARY 


Land managers often use plot data to assign pine bark beetle hazard ratings for lodgepole pine (LPP), 
ponderosa pine (PP), western white pine (WWP), whitebark and limber pines (WBkP-LmP), and all pine 
species (ALLPINE) (Randall et al., 2011). The USDA Forest Service Northern Region (R1), developed 
pine bark beetle hazard ratings (very low, low, moderate, or high) for lodgepole pine (LPP), ponderosa 
pine (PP), western white pine (WWP), whitebark and limber pines (WBkP-LmP), and all pine species 
(ALLPINE) (Randall et al., 2011). These hazard ratings are available for all inventory data, both stand 
exam and Forest Inventory and Analysis (FIA) plot data residing in the National Forest Systems 
Vegetation Database, FS Veg and for growth and yield simulations using the Forest Vegetation Simulator 
(FVS). 


Since stand exam data can be attributed with pine beetle hazard, this information can be displayed with 
the associated stand polygon layer. However, because many stands have not been examined, or stand 
exam data is outdated, hazard maps may not be contiguous or current. It was desired to produce a spatial 
layer of pine beetle hazard rating, for the Region that could be used by land managers. 


The Northern Region has an existing vegetation classification system (Barber et al., 2011) which 
determines dominant tree species, tree size class, and canopy cover class. This classification system 
meets national tech guide standards as defined in the USDA Forest Service Existing Vegetation 
Classification, Mapping, and Inventory technical guide (Citation, the one below is old). This 
classification system be applied to inventory data and is used to derive the existing vegetation spatial 
database for the Northern Region. The Northern Region Vegetation Mapping Project (R1-VMap) 
(Brewer et al., 2004) is the most current, region-wide spatial database depicting existing vegetation. This 
spatial dataset displays the following vegetation attributes: species dominance type, tree size class, and 
tree canopy cover class. Rl1-VMap meets, and in many cases exceeds, the requirements of an existing 
vegetation system as defined in the USDA Forest Service Existing Vegetation and Mapping Technical 
Guide (Nelson et al., 2015) and the National Vegetation Classification Standard (FGDC NVC, 2008). 
R1-VMap includes lands of other ownerships near and within National Forest boundaries, and can be 
used to create specialized map products useful to land managers 


A spatially balanced sample of inventory plots (n=3324) were used to explore the relationship of pine 
beetle hazard to R1-ExVeg classifications for dominance group, size class, and canopy cover class with 
the goal of assigning pine beetle hazard to Rl1-VMap. Based on our analysis we determined that R1- 
ExVeg attributes of dominance group 60% plurality (Dom60) and tree size best predicted pine bark beetle 
hazard rating class. Using a standard rule base we assigned the beetle hazard rating that best fit each 
Dom60 / tree size class in the RI-VMap which provides a region-wide pine beetle hazard spatial 
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depiction. In addition, we explore options for customizing this process for individual National Forests. 
Other methods could have been used to assign hazard to plots and/or R1-VMap. 


The resulting beetle hazard maps are subject to the same strengths and weaknesses inherent in the base 
R1-VMap data. Thus, care should be taken in using them at appropriate spatial scales. In addition, 
hazard ratings, in themselves, do not provide a direct measure of tree mortality, although it may be 
possible in the future to predict tree loss in RI-VMap polygons, based on their hazard rating, particularly 
for mountain pine beetle at epidemic levels. 


INTRODUCTION 


Pine bark beetles, predominantly mountain pine beetle (MPB) (Dendroctonus ponderosae Hopkins) and 
to some degree western pine beetle (WPB) (Dendroctonus brevicomis LeConte), are the most destructive 
bark beetles (Coleoptera: Scolytidae) in western North America and may conflict with pine species 
management objectives. For significant pine bark beetle-caused mortality to occur, three things are 
necessary: 1) suitable host trees, 2) pine bark beetle populations, and 3) favorable weather conditions. Of 
these three the only factor land managers can influence is the amount and arrangement of suitable host 
trees. A landscape-level map identifying forested areas most susceptible to significant pine bark beetle 
impacts can help Northern Region (R1) National Forest managers determine where limited management 
resources can be dedicated. 


Pine bark beetle hazard rating systems which use stand characteristics such as basal area (BA), average 
tree diameter at breast height (DBH), and percent of susceptible host species in the stand, enable land 
managers to identify forested areas at high risk to significant pine-volume loss where current stand 
information is available (see Randall et al. 2011). As stand exams are conducted less frequently, the area 
covered by and the accuracy of stand-based hazard rating maps has diminished. 


The Northern Region Vegetation Mapping Project (R1-VMap) (Brewer et al., 2004) is the most current, 
region-wide existing vegetation spatial layer. This spatial dataset contains vegetation attributes as defined 
by the Northern Region’s Existing Vegetation Classification System (R1-Ex Veg) (Barber, Bush, & 
Berglund, 2011). R1-VMap meets, and in many cases exceeds, the requirements of an existing vegetation 
system as defined in the USDA Forest Service Existing Vegetation and Mapping Technical Guide 
(Nelson et al., 2015) and the National Vegetation Classification Standard (FGDC NVC, 2008). R1-VMap 
includes lands of other ownerships near and within National Forest boundaries, and can be used to create 
specialized map products useful to land managers. Our goal was to take R1-VMap and create a protocol 
that allows managers to assign pine bark beetle hazard, useful in management decisions. 


METHODS and RESULTS 
DATA RESOURCES 
R1 Existing Vegetation Classification System (R1-ExVeg) and Vegetation Map (R1-VMa 


The Northern Region Existing Vegetation Mapping Program (VMap) database and map products provides 
the Northern Region with a geospatial database of existing vegetation produced using consistent 
analytical methodology according to the Existing Vegetation Classification and Mapping Technical Guide 
(Nelson, et. Al. 2015) to support the Region 1 Multi-level Classification, Mapping, Inventory, and 
Analysis System, R1-CMIA (Berglund et. al., 2009). The VMap database provides four primary map 
products; lifeform, tree canopy cover class, tree size class, and tree dominance type to support mid- and 
base-level analysis and planning. VMap uses the Region | Existing Vegetation Classification System 
(R1-ExVeg) (Barber, et.al. 2009) in its map unit design. The system was designed to allow consistent 
applications between Regional inventory and map products within the R1-CMIA framework. VMap is a 
remote sensing derived product. As such, it uses a combination of high resolution airborne imagery and a 
nationally available digital elevation model (DEM). An accuracy assessment is conducted to provide a 
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validation of the data, giving an indication of reliability of the map products, so that managers are fully 
informed throughout the decision making process. Estimates of overall map accuracy and confidence of 
individual map classes can be inferred from the accuracy assessment error matrix derived from the 
comparison of known reference sites to mapped data. 














Ri Summary FIA Plot Inventory Database (R1-SDB) 


Forest Inventory and Analyses (FIA) plot data (USDA Forest Service 2016) provides a spatially balanced, 
statistical sample appropriate to use for broad level planning and analysis (Bush, 2014). It covers all 
lands regardless of land-use designation or management history, making it representative of vegetation 
across the Region. The FIA inventory used for this hazard rating was reasonably current with plots 
installed by FIA between the 1993 and 2002 period (periodic inventory plots). In 2003, Interior West 
(IW) FIA started an annual inventory which re-measures 10% of the plots every year. These newer plots 
were not available during model development. The periodic data in the Northern Region’s FS-VEG RI- 
SDB database (R1-SDB ver. 2003) covered 3423 forested FIA plots located on Forest Service lands. 
However, 99 plots had no vegetation type (all were “non-forest”) and were dropped so that 3324 plots 
were used in our analyses. 


Pine Bark Beetle Hazard Ratings Using Plot Data (Plot Hazard) 


Pine bark beetle hazard ratings for five different pine hosts were assigned to the FIA plots. The criteria 
for assigning each of the pine bark beetle hazard rating classes (Table 1) to inventory data is summarized 
in Randall and others (2011). Each plot’s pine bark beetle hazard ratings were accurate as of the time of 
plot inventory. All plots were surveyed prior to the current MPB outbreak, and hazard ratings only 
needed to apply to the plot for the date data was collected, no attempts were made to “grow” plots from 
the date of the last inventory to the present. 


Abbreviations for the five pine host models are: 
e ALL PINE: Pine bark beetles in all pine species 
e LPP: Mountain pine beetle in lodgepole pine 


e PP: Bark beetles in ponderosa pine (includes mountain pine beetle, possible Jps species, and 
western pine beetles west of the Continental Divide). 


e WWP: Mountain pine beetle in western white pine 


e WBkP-LmP: Mountain pine beetle in whitebark and limber pines 






_ DESCRIPTION 
| EOS ah ~ VERY LOW HAZARD: host trees are not present in the plot data; pine bark beetle impact is 
expected to be minimal or non-existent 
: 1 piace. LOW HAZARD: - host trees are either not abundant and/ or of sufficient < size eto be 
susceptible to pine bark beetle attack resulting in pine mortality 
eal MODERATE F HAZARD: plot data indicates that there are enough host trees of the 
| appropriate size to be susceptible to pine bark beetle attack and subsequent pine 
mortality 
aa. HIGH HAZARD: “plot data indicates that there are enough host trees of the appropriate 
‘size to be susceptible to pine bark beetle attack and d significant subsequent pine ‘mortality ‘| 


4 NO HAZARD: non- forested lands 


The total number of FIA plots in each of the four hazard classes for each of the five pine host models is 
given in Table 2. 


Table 2. Number of plots in each hazard class by pine host hazard rating 


i oo | 
eal a line |e 
acl ee ee 
CTorAL 

















885 
521 SAL 
256 
824 44 
TOTAL 


*Western white pine is modelled for west of Continental Divide only 





Vegetation Attributes (R1-ExVeg) 


Each FIA plot was classified into dominance type, tree size, and tree canopy cover classes using the R1- 
ExVeg classification system (Barber, Bush, & Berglund, 2011). These classifications are described 
below. 


Dominance Type 
R1-ExVeg dominance type describes tree species composition in a forest. 


e Dominance group 6040 (DG6040) is based on whether the dominant tree species constitutes 
>60%, 59-40%, or <40% of the stand (based on logic using one or more factors of canopy 
cover, BA, or trees per acre but for stands with >=20 ft? total BA dominance is usually based 
on relative BA). While DG6040 is the finest thematic dominance group classification system 
mapped in R1-VMap, mapping personnel did not feel there was an adequate accuracy 
assessment for this classification system for mid- and broad-level assessments. For this reason 
we did not consider DG6040 as the vegetation dominance group classification for the R1- 
VMap hazard rating system. 


e Dominance 60% plurality (Dom60) classes include only single-species classes and mixed- 
species classes. This creates a map or inventory compilation with classes that are based on > 
60% abundance of an individual species and three heterogeneous mixed species classes. 
Dom60 classes are most useful when the management question of interest requires relatively 
pure vegetation types. 


e Dominance 40% plurality (Dom40) classes consolidate all single species classes and single 
species-mixed species classes together based on the dominant species present. This creates a 
map or inventory compilation with classes that are based on >40% abundance. 


Size Class 


Tree size is a classification of the predominant DBH of live trees within a plot. Size class for 
inventory data is based on basal area weighted average diameter. Each plot is assigned one of four 
tree size classes: 


e Seedling- <4.9” DBH 
5.0”-9.9” DBH 
10.0”-14.9” DBH 
>15”+ DBH 




















Canopy Cover 


Tree canopy cover describes the proportion of the forest floor covered by the vertical projection of the 
tree crown. . For inventory data, canopy cover is calculated using the Forest Vegetation Simulator 
(FVS)Each plot is assigned one of five tree canopy cover classes: 


0-9.9% 
10-24.9% 
25-39.9% 
40-59.9% 
60% + 


Table 3. National Forests (and 
identification number) in the 

Northern Region identified as 
principally east or west of the 
Continental Divide 


EAST SIDE WEST SIDE 


Beaverhead- 


Lewis & Clark (15) | Lolo (16) 
enareeeeited|Nez Perce'(17) 







Continental Divide 


Past and current climate, vegetation, and insect activity can differ 
based on whether you are on the east or west side of the 
Continental Divide (CD). In general, Forests in the Northern 
Region reside on one side or another of the CD. For a few 
Forests some boundaries do cross over this line, however, this is 
often minimal. Thus, plots were assigned as east or west based 
on which Forest they occurred on (Table 3). 








PROCESS 


In existing hazard rating systems, hazard rating classes are assigned to stands based on plot level data. 
R1-VMap polygons do not have plot level data, but they do have R1-ExVeg classes that can be used as 
surrogate stand descriptions. In order to create pine bark beetle hazard rating maps using R1-VMap, pine 
beetle hazard ratings had to be assigned to one or a combination of R1-ExVeg class attributes available in 
R1-VMap. 


Determining R1-ExVeg Attributes Best for Modelling 


The first step in creating pine bark beetle hazard rating maps using R1-VMap required analyses of FIA 
plot inventory data (R1-SDB) to determine which R1-ExVeg vegetation attributes best predict beetle 
hazard rating. 


Through graphing and summarizing exercises we found moderate and high pine bark beetle hazard ratings 
were concentrated in fewer Dom60 classes than Dom40 classes, although there were an equal number of 
classes in DG60 and DG40 represented in the FIA plot data. Overall results indicated that the DG60 
classification system did a better job of identifying stands with significant host components than DG40. 


We then used the R statistics program (R Development Core Team, 2011), rpart module (Therneau and 
Atkinson, 2011) to evaluate R1-ExVeg variables of Dom60, size class, canopy cover, and CD for their ~ 
importance in predicting pine bark beetle hazard rating classes. Results of rpart showed Dom60 to be the 
primary variable of importance in predicting a plots pine bark beetle hazard class for all four individual 
host models, followed in importance by tree size class. If tree canopy cover was determined to contribute 
to explaining hazard, it was usually the 3" factor in the host model (Table 4). 


Continental Divide was identified as the primary variable for the WBkP-LmP model (Table 4). However, 
it is likely the importance of Continental Divide is an artifact of MPB-caused mortality in this host type 


if 


west of the Divide several decades prior to collection of FIA periodic data. We do provide both an East 
and a West version of the WBkP-LmP hazard model (Appendix C) that Forests may explore to determine 
if these better depict what is expected in their area. The Regional model, however, is recommended for 
use if close comparison of models is not done. 





Table 4. Order of importance of variables used to build the best model for 
predicting pine beetle hazard for 4 habitat grouping. Variables given for the first 3 
splits, all branches. An ‘X’ indicates the branch dead ended; ‘n/a’ indicates no 
applicable data available (branch terminated previously). 


HAZ RATING 1° SPLIT 2"¢ SPLIT 3"? SPLIT 


Size Class x 
ALL PINES Dominance Type i eae ees Size Class. 
Size Class x 
Dominance Type 







































X (n/a) 

LODGEPOLE Dominance Type {ofa} 
Dominance Type Size Class 
Size Class 






(n/a) 
(n/a) 

X 
Canopy Cover 
X (n/a) 
Continental Divide (nfa) 
Dominance Type X 
Dominance Type 









Dominance Type 






- PONDEROSA Size Class 
















WHITEBARK / LIMBER 






Assessing Model Error Rates 
Overall Model Error 


Regional rating of plot-level data for pine bark beetle hazard is based on considerable published literature 
(See Randall et al., 2011) and makes use of a range of detailed inventory data. However, data currently 
available in R1-VMap is much less detailed, with categorized variables that may not reflect divisions that 
are biologically relevant to pine bark beetles. Generalization of the rating system to a few categorical 
factors, when hazard is based on more detailed stand information should be recognized as a significant 
source of error. Thus, it is be important to understand how well these R1-VMap values serve as 
substitutes in describing plot-level hazard. 


In addition to using the R statistical package, RandomForest (Liaw and Wiener, 2002) was used to further 
test whether the two-factor model we had chosen using rpart was the most appropriate model using 
overall model error rates. Comparison of possible models using R module RandomForest supports the 
importance of Dom60, followed in importance by Size Class and Canopy Cover (Table 5). Across all five 
host models and for ALL PINES in particular, the best performing variable model with least number of 
factors was usually Dom60/Size Class. In a few cases other models performed equally well. 





Table 5: Error rates (rounded to closest 1%) for all possible variable models for each’ 
pine host hazard rating. Values in blue print indicate models with lowest errors 
| or equal-to-lowest error with fewest variables (simplest model). 


VARIABLE ALL — Lodgepole Ponderosa Whitebark Western 
MODEL PINE / Limber White 
Dom60 49% 38% 11% 23% 15% 
| Dom60 + 46% 35% 10% 23% 15% 
| Canopy 
: Canopy + Size 55% 45% 13% 25% 15% 
Dom60 + Size 43% 33% 10% 23% 15% 
: Dom60 + Size + 43% 34% 10% 23% 15% 
| Canopy 
| Dom60+Sizet+ 43% 34% 11% 24% n/a 
Continental 


Divide 











Table 6: Individual hazard rating error 
rates for all pine species (ALL PINE) 
for each of the six variable models 
tested. See text for more detail on how 
to interpret numbers. 


ALL PINES 


PREDICTED HAZARD 
DomGroup 60% Pluralit 


274 0 32% 

21 179 104 96% 

8 528 185 52% 
7 252 


LoS 23 66% 
30 498 203 54% 


203 29% 


302 16 39% 
192 eee OL OS 59 79% 
26. 2998 is VPS AST 


188 172 
34 604 
eeoes 404 
DomGroup 60% Plurality*Size Class 
*Canopy Cover 


ACTUAL PLOT HAZARD RATING 


37% 


ity/al 20 70% 


587 4.173 46% 

Me 34 0210 580 30% 

DomGroup 60% Plurality*Size Class 
*CanopyCover*Continen 
51 227 

69% 

51% 











Error by Hazard Class: 


RandomForest was also used to provide prediction error 
rates for individual hazard rating classes, i.e. how well did 
our models using R1-ExVeg categories predict the plot 
hazard ratings. So, whereas the ALL PINE hosts model 
using Dom60/Size Class as factors had an overall error of 
43% (Table 5), individual hazard rating level errors are 
28%, 66%, 54%, and 29% for hazards of VERY LOW (0), 
LOW (1), MODERATE (2), and HIGH (3), respectively 
(Table 6). 


To understand the confusion matrix reported by 
RandomForest output (Table 6), consider that each row 
represents all FIA plots within the hazard rating class given 
on the left (based on plot variables). Along the row, 
individual columns show how many of those plots the new 
model (using vegetation classifications for Dom60 and Size 
Class) predicted would be in each hazard rating. For 
example, in the model shown in Table 6 using only Dom60 
to predict hazard, of the 824 FIA plots with a HIGH (3) 
hazard rating (33+7+252+532; also see Table 2), the 1- 
variable model correctly predicted 532 plots as being High 
hazard, placing most of the miss-categorized plots 
(33+7+252) into MODERATE (2) hazard (252). 


Review of the individual hazard level ratings suggests that 
our preferred two-factor model of Dom60/Size Class does 


- better at predicting areas with VERY LOW and HIGH 


hazards, which are usually the most important categories 
for land managers (Appendix A). LOW and MODERATE 
categories have higher error rates and should be used with 
caution. 
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Assigning Pine Bark Beetle Hazard Ratings to R1-VMap Polygons 


Using the two-variable model that best predicted plot hazard for all five beetle host types, we then had to 
determine which hazard rating to assign each Dom60/Size Class vegetation group. Our process for 
assigning pine bark beetle hazard ratings (Table 1) to each combination of Dom60/Size Class used the 
logic outlined in Table 7. This logic is based on professional judgment, and efforts to reflect hazard 
ratings of ground plots that might have had similar dominance type and tree size class. 





1) If the percent of inventory plots with Very Low hazard (0) is >75%, then Dom60/Size Class hazard 
__is Very Low (0); 
| Else go to 2 


2) ‘Tf the percent it of in inventory plots with. High hazard (3) i is >50% and >80% are in High (3) OR 
Moderate (2), then Dom60/Size Class hazard is High @)s_ 
| Else- go to 3 


| 3) If the percent of ir inventory plots with High(3) OR Moderate (2) hazard is >40%, then Dom60/Size 
Class hazard is Moderate (2); 


Else Dom60/Size Cl: Class ha: hazard is Low v (1) 
_(Non- forested areas é areas are giv given a hazard rd category of 4) 


wee 


While we had a large number of FIA plots for some Dom60/Size Class combinations, some combinations 
had few or no plots and we were forced to use expert opinion to assign a hazard rating. We wanted to 
capture the sample strength for the pine bark beetle hazard ratings assigned to each Dom60/Size Class 
combination. Thus, we also provide a sample strength indicator based on the number of inventory plots 
within each vegetation grouping as outlined in Table 8. These codes are used throughout Appendices B 
and C and are indicated by the color of the data cell. 


Table 8: Description of Sample Strength Indicators (and color codes used in appendices) 

1 GOOD SAMPLE SIZE: FIA sample was sufficient to assess the hazard for vegetation o 
group (n>10 FIA plots in vegetation group). (no color) . aoxer 

2 SMALL SAMPLE SIZE: FIA sample was not sufficient to assess the hazard for vegetation 
group (n_<10 FIA plots in vegetation group); Expert call supported by sample. (gold) _ 

3 SMALL SAMPLE SIZE: FIA sample was not sufficient to assess the hazard vegetation 
group (n<10 FIA plots in vegetation group); Expert call different from sample. (gold) 

4 NO SAMPLE: No FIA sample plots in DG60/ Size Class combination; Expert call on 
hazard. (orange) ant’ 

5 NON-FORESTED: vegetation group was non-forested, so no hazard assigned. (N/A) 











Attaching the hazard rating and sample strength codes to the R1-VMap polygons allows spatial depiction 
and analyses. Figure 1 shows what the ALL PINES hazard ratings are for Region 1, with a close up that 
also shows how sample strength has been presented in the current map available from Region | 
Engineering, Geospatial staff. 
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DISCUSSION 


EXPLORING ERROR OF HAZARD ASSIGNMENTS 


The database of FIA plots with plot level hazard ratings was updated to include the new hazard rating 
corresponding to the Dom60/Size Class of that plot. With this, we further explored sources of error and 
error rates with the goal of possible refinement of the hazard rating model. 


Comparing R1-VMap and R1-SDB for Forest-level hazard 


We conducted an additional evaluation of the 3324 FIA plots to compare the estimates of hazard using 
Randall and others (2011) to the estimates of hazard using the equivalent R1-Veg classifications. .Table 9 
shows how well the two methods of hazard rating correspond. Note VERY LOW (0) and LOW (1) 
hazards were combined to smooth out the high variability between datasets for these two categories. 

*Data suggests VERY LOW (0) and LOW (1) hazard ratings should be combined to provide most similar area estimates 


Table 9: Comparison of area predicted to be in each of three hazard rating classes based on hazard 
ratings assigned using either Randall and others (2011) or R1-Veg Dom60/Size Class. The ALL PINE 
pcmeerrroes results were used for all 3324 plots.. 


BASED ON RANDALL ET AL BASED ON DOM60/SIZE CLASS DIFFERENCES BETWEEN THE TWO 
ee RATINGS laa RATING ANE 
Forest Name bse 2 i ey: Hie 
(Number) Plots low re i 

















soa [| 
io |= ae 
Paieevoorts) | 220 | 
Prankanaee) | 369 | 
Feearwaer(5) | 205 | 
Fraheastao) | 295 | 
rooteai (gy | 344 | 
Fisiowy | a2 | 
[nex Pere(a7h | 25 | 
3324 | 


**In most cases, estimates differed iby 10% c or ress TcepE Ere noted i in ied 















































Across the 12 Forests and Regionally the number of plots in VERY LOW+LOW, MODERATE, and 
HIGH hazard for both hazard rating methods are within 10% agreement. Five Forests have greater than 
10% disparity in the MODERATE rating, with one Forest (Kootenai) also off >10% in the lowest hazard 
rating. In no case was disparity greater than 15%. Overall, this is a good level of agreement at the 
landscape level. 


Note that this comparison was done on the older inventory plots. The depiction of error between the two 

hazard rating systems (Randall et al 2011 versus Dom60/Size Class) is appropriate, but the numbers are 

not indicative of current Forest-wide conditions. To obtain more current Forest-wide conditions, more 

recent inventory data (preferably after bark beetle, fire or other large disturbances were completed) would 

have to be used. Use of our rating system on newer R1-VMap coverages is discussed in the section on 
Exploring Error: Age of Data. 
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Note that VERY LOW and LOW hazard ratings were combined to obtain similar prediction levels. Stand 
level plot hazard ratings placed approximately 2/3 of VERY LOW and LOW hazard plot into VERY 
LOW with the other 1/3 into LOW (Regional estimate; Forest estimates averaged 3/5 and 2/5, 
respectively). Distribution of these lower hazard ratings based on R1-Veg classification put most plots 
into the LOW hazard with approximately 1/20 and 19/20 assigned to VERY LOW and LOW 
(respectively) both for the Region and as the Forest average. This may suggest that in our decision matrix 
(Table 6) our requirement that, for any Dom60/Size Class vegetation group, >75% of plots have a rating 
of VERY LOW (0) for the vegetation group to be classified as having VERY LOW (0) hazard is overly 
stringent; perhaps >51% would be more appropriate. 


Classification error at plot level (mismatch error) 


If we consider the hazard rating assigned to the RI-SDB FIA plots using Randal and others (2011) as 
being the best estimate of “real” MPB hazard and the hazard assigned to the R1-ExVeg classification 
(used in R1-VMap) as “estimated” or modelled hazard, a comparison of the two ratings will highlight 
Dom60 groups where errors are highest. Table 10 provides a summary view of mismatch in these hazard 
ratings. Of greatest concern are 1) errors in common Dom60 vegetation types and 2) errors greater than 
one rating of mismatch. Also, depending on the intended use of the map, modelled ratings that are lower 
than what might be expected on the ground may be of greater concern than if the model over estimates 
hazard in some polygons or plots. 


Where the Forest-wide evaluation in Table 9 allowed higher and lower ratings to cancel each other out 
over the landscape, this evaluation of error is less forgiving. In general, we found that our modelled 
hazard ratings provide a correct match to the plot hazard rating 46% or about half of the time. 
Approximately 36% or one third of the time the model predicted a higher hazard than the FIA plots were 
given, and 18% of the time the model gave a lower hazard than the FIA plots. Taking the conservative 
approach and accepting where the model slightly overestimates hazard (i.e. mismatches of -1) we looked 
at the two situations that could cause the greatest concern to managers. 


1. Mismatches where the model overestimated hazard by more than one level (i.e. mismatch of -2) 
that occurred more than 10% of the time (e.g. PIEN, PIFL2, PSME, TMIX in red highlight), and 

2. Mismatches where the model underestimated hazard by any amount more than 25% of the time 
(e.g. ABGR, ABLA, IMIX in yellow highlight). 


If we further explore only Dom60 types that make up a significant portion (>10%) of the landscape (in 
this case, the Region), we are left with IMIX, PSME, and TMIX as the vegetation types with the most 
errors of concern over the largest area. These three vegetation types make up approximately 56% of the 
FIA plots for the Region. Dom60 types with fewer than 10 plots do not have sufficient information for 
further dissection and make up too little of the landscape to evaluate at this level (grey highlight). 


Where the Forest-level evaluation in Table 9 shows which Forests may benefit from a reevaluation of the 
models, this evaluation of mismatching hazard ratings by Dom60 suggests which vegetation types warrant 
closer inspection to determine if additional factors might be considered to make more accurate hazard 
ratings. Note that a significant portion of single level model errors (-1 and 1) may be between VERY 
LOW (0) and LOW (1) hazard levels which have been combined in Table 9. 
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Table 10. Mismatch between R1-SDB FIA plots (“reality”) and R1- 
Veg category rating assigned to R1-VMap (“modelled”) using the 
ALLPINE model. 


modelled as higher | —*||__realityhigher hazard | Number of 
hazard than reality* than modelled** 


ale Aa 


aN iematchas' where the model. overestimated hazard by more than one level (i.e. 
mismatch of -2) that occurred more than 10% of the time are highlighted in red 
**Mismatches where the model underestimated hazard by any amount more than 25% of 
the time are highlighted in yellow 





Bimodal data 


One important source of mismatch error is a binomial distribution of the plot hazard ratings. For 
example, several Dom60 cover types had plot hazard ratings of 0 or 2 (Fig. 2), but when a hazard rating 
was assigned to the Dom60 group the decision tree resulted in assignment of hazard rating 1. Although 
this approximates an “‘average” of plot ratings, it results in a high rate of hazard misclassification. 


Regionally, this misclassification was particularly strong in PSME and TMIX. We conducted a cursory 
evaluation of several additional factors at the Regional level that might help divide stands with VERY 
LOW (0) from MODERATE (2) hazard. These factors included elevation, aspect, slope and CD. 
Identification of a factor that predicted placement of plots into one or the other hazard rating would 
greatly decrease misclassification errors. Unfortunately, none of the most likely factors evaluated 
appeared to be useful. Assessment was visual (graphics) and a full statistical analyses of potential 
variables may. be needed. However, we feel this would be most appropriate at the Forest level if 
warranted. 


1) 








mm VERY LOW (0) 
m LOW (1) 

- MODERATE (2) . 
m= HIGH (3) 











IMIX PICO 


n 
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xe) 
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eo 
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Figure 2. All Dom60 vegetation types for the Region with >100 FIA plots in the R1- 
SDB. Several types show a bimodal distribution of 1/3 or more of plots in the VERY 
LOW (0) and MODERATE (2) hazard classes (circled in red). PSME and TMIX are of 
greatest impact to overall model error rates due to the high number of plots in these 
groups. This data is for the ALL PINE model. 
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Age of Data 


Since our first models were completed the Region updated the R1-SDB dataset of FIA plots. This new 
dataset (ver. 2007) has approximately half of the plot data replaced by information gathered since 2003. 
Using the same Dom60/Size Class model, and the thresholds used in our original calculations (Table 7) 
we recalculated the ALLPINE, LPP, and PP models. The updated calculations resulted in some changes 
in some hazard ratings. Many of these were between VERY LOW (0) and LOW (1) which was 
considered of little importance, especially since these two ratings might best be combined for managers. 
However, there were a couple of changes that did appear important. The updated LPP model assigned 
PICO 5-9.9° and PICO 10-14.9” to MODERATE (2) hazard rather than HIGH (3). This same adjustment 
was suggested for PIFL2 5-9.9” in the ALLPINE model. 


These suggested hazard rating shifts likely reflect recent changes in stand conditions due to the mountain 
pine beetle (MBP) outbreak. Because Canopy Cover was indicated by statistics to be a relatively 
unimportant factor in the hazard rating model developed using the largely pre-outbreak data, the recent 
drop in stand density from MPB-caused mortality is not reflected in the Dom60/Size Class model. 
Although average tree diameter will decrease as MPB remove the larger trees, many LPP stands are 
similar enough in size that most stems fit within one of the Size Class categories, especially the 5-9.9” 
DBH Size Class. Thus, plot-level hazard rating that uses stand density (Randall et al., 2011) may now be 
lower after MPB activity, yet the hazard rating of that same plot using only Dom60 and Size Class 
assignment would stay the same. . 


Ultimately, this means that R1-VMap coverages created after significant MPB-caused mortality has 
occurred will likely show high hazard remaining in stands that have already sustained MPB activity. As 
with use of older R1-VMap coverages (pre-MPB outbreak) this stand density alteration due to MPB 
activity would need to be accounted for to represent current conditions. With new R1-VMap products 
coming on board we will have to re-create our hazard ratings which may require use of a model that 
include Canopy Cover or a completely different approach all together. 


Updating R1-VMap 


VMap is a geospatial database that allows for editing of polygon labels when better local information 
exists to do so. Updates and corrections to the R1-VMap product will increase accuracy of landscape pine 
bark beetle hazard depiction. Forest or analyses area maps that have VMap polygons edited to account 
for better local data sources can have the Regional (or updated) hazard matrix attached. In fact, our 
hazard matrix can be applied to any spatial dataset that uses R1-Ex Veg classifications of Dom60 and tree 
size class. 


INDIVIDUALIZING HAZARD RATINGS 


It is important to note that a more complicated model of hazard rating for use in Rl-VMap is not indicated 
by our evaluation of the data. Much of the error in the models is due to the R1-ExVeg classification 
system using categories that do not reflect environmental thresholds biologically important to MPB. Until 
more appropriate divisions can be obtained for modeling, a high level of misclassification of bark beetle 
hazard in RI-VMAP will remain. - 


However, a model more specific to a given area or Forest (versus Region) may be obtained by 
recalculating hazard ratings using only plots from the Forest. There are also some predominant Dom60 
cover types that would benefit from consideration of additional factors (e.g. Canopy Cover) or that have a 
binomial distribution where one of the two predominant hazard ratings would be more appropriate than 
the middle or average rating. 


Two examples of re-evaluations are given in Appendix D. The Gallatin NF example looks at Forest- 
specific alterations to the Regional hazard ratings. This was prompted by the 15% error in the 
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MODERATE hazard type across the forest (Table 9). The Lolo NF example looks at the R1-VMap 
polygons within an analyses area to explore uncertainty in the hazard rating of the IMIX cover type. 


Although not explored in any example, our evaluation of the updated FIA data suggests that for PIFL2 5- 
9.9” in the ALLPINE model, and PICO 5-9.9’ and PICO 10-14.9” in the LPP model the HIGH (3) hazard 
rating may need to be degraded to MODERATE (2) to reflect recent MPB activity that would have 
lowered stand density below the high-hazard threshold. One way to address this would be to consider 
Canopy Cover in these categories. For example, if canopy cover is <60% (or <40%) HIGH hazard is 
maintained, but if below that threshold the hazard is dropped to MODERATE. 


If Forests want to change ratings they should work with FHP and be prepared to document reasons for all 
changes to maintain repeatability and supportability should results be questioned. 
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APPENDIX B: Hazard ratings and plot counts for five host models using Dominance Group 
60% Plurality / Size Class categories 


BARK BEETLES, PRINCIPALLY MOUNTAIN PINE BEETLE, 
In ALL PINE HOST SPECIES (ALL PINE) 







































































































































































































HAZARD RATING* i SAMPLE SIZE* 
Dominance Size Class (inches DBH) Dominance Size Class (inches DBH 
Group 60% | 0.0-4.9 05.0-09.9 | 10.0-14.9 | 15.0+ Group 60% | 0.0-4.9 | 05.0-09.9 | 10.0-14.9 | 15.0+ 
Plurality! Total Total Total Plurality' Total Total Total Total 
ABLA 97 100 39 
BEPA 
CELE3 
IMIX 
JUNIP 
LALY 
LAOC 
PIAL 30 
PICO Poe 
PIEN 13 
PIFL2 
PIMO3 
PIPO PIPO 
POPUL POPRBE 
POTRS POTRS 
PSME PSME 
THEE 3] 
TMIX 185 
TSHE 1] 
| TSME 13 











*Color coding was used to indicate sample sizes; estimates based on >10 plots has no color; on 1-9 plots are gold, and on 0 plots are orange. 
Hazard ratings based on <10 plots were determined using expert opinion. For calls based on 1-9 plots, changes to predicted hazard based on 
professional opinion are indicated as the suggested hazard with predicted in () following [e.g. “1 (not 0)”]. 


'Dominance Group 60% Pluralities of HMIX, FRPE, and TABR2 had no plots available for analyses and are not expected to have sufficient 
MPB host to warrant a hazard above 0 so are not included in this table. 


APPENDIX B (cont.) 




















MOUNTAIN PINE BEETLE 
In LODGEPOLE PINE (LPP) 
HAZARD RATING* SAMPLE SIZE* 
Dominance Size Class (inches DBH) aa Dominance Size Class (inches DBH 
Group 60% | 0.0-4.9 | 05.0-09.9 | 10.0-14.9 | 15.0+ Group 60% 05.0-09.9 | 10.0-14.9 
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Total 


: Total Total 















































































































































*Color coding was used to indicate sample sizes; estimates based on 

>10 plots has no color; on 1-9 plots are gold, and on 0 plots are orange. Hazard ratings based on <10 plots were determined using expert 
opinion. For calls based on 1-9 plots, changes to predicted hazard based on professional opinion are indicated as the suggested hazard with 
predicted in () following [e.g. “1 (not 0)”]. 


'Dominance Group 60% Pluralities of HMIX, FRPE, and TABR2 had no plots available for analyses and are not expected to have sufficient 
MPB host to warrant a hazard above 0 so are not included in this table 
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APPENDIX B (cont.) 





















































































































































































































































PINE BARK BEETLES* 
In PONDEROSA PINE (PP) 
HAZARD RATING** | SAMPLE SIZE** 
Dominance Size Class (inches DBH Dominance Size Class (inches DBH) 
Group 60% | 0.0-4.9 | 05.0-09.9 | 10.0-14.9 | 15.0+ Group 60% | 0.0-4.9 | 05.0-09.9 | 10.0-14.9 | 15.0+ 
Plurality! Total Total Total Total Plurality! Total Total Total Total 
ABGR ABGR 10 34 65 
ABLA we ABLA 97 100 39 
BEPA BEPA 
CELE3 CELE3 
IMIX IMIX 
_JUNIP JUNIP 
LALY LALY 
LAOC LAOC 10 10 
PIAL PIAL 30 43 16 
PICO PICO re. 808 572 3s se 
PIEN PIEN 13 44 5] 
PIFL2 PIFL2 
PIMO3 PIMO3 
PIPO PIPO 41 
POPUL POPUL 
POTRS POTRS 
PSME eaieep on enw 00 | Po!" ||PSME 229 
THPL 0 THPL 31 
TMIX roy rater 0.) Eto | 0} | EMIX 
| TSHE 0 TSHE elit 
0 TSME 13 





*Includes activity by mountain pine beetle, western pine beetle (only present west of Continental Divide) and /ps beetles 

**Color coding was used to indicate sample sizes; estimates based on >10 plots has no color; on 1-9 plots are gold, and on 0 plots are orange. 
Hazard ratings based on <10 plots were determined using expert opinion. For calls based on 1-9 plots, changes to predicted hazard based on 
professional opinion are indicated as the suggested hazard with predicted in () following [e.g. “1 (not 0)”). 


'Dominance Group 60% Pluralities of HMIX, FRPE, and TABR2 had no plots available for analyses and are not expected to have sufficient 


MPB host to warrant a hazard above 0 so are not included in this table. 


APPENDIX B (cont.) 
MOUNTAIN PINE BEETLE 
In WHITE BARK and LIMBER PINES (WBkLmP)* 





HAZARD RATING** SAMPLE SIZE** 
















Dominance Size Class (inches DBH Dominance Size Class (inches DBH 


















Group Group 
60% 0.0-4.9 05.0-09.9 | 10.0-14.9 15207, 60% 05.0-09.9 | 10.0-14.9 | 15.0+ 
Plurality! Total Total Total Total Plurality! Total Total Total 





































































































































































ABLA 97 100 

BEPA 

CELE3 

IMIX 
JUNIP 

LALY 

LAOC 

PIAL PIAL 30 43 

PICO we uameon 372 138 

PIEN PIEN 13 44 

PIFL2 PIFL2 

PIMO3 PIMO3 

PIPO 
POPUL POPUL 

POTRS POTRS 

PSME PSME 
THPL THPL 
TMIX 
TSHE 
TSME 0 





*See Appendix C for details on possible model alterations for east and west of the Continental Divide 


**Color coding was used to indicate sample sizes; estimates based on >10 plots has no color; on 1-9 plots are gold, and on 0 plots are orange. 
Hazard ratings based on <10 plots were determined using expert opinion. For calls based on 1-9 plots, changes to predicted hazard based on 
professional opinion are indicated as the suggested hazard with predicted in () following [e.g. “1 (not 0)”). 


'Dominance Group 60% Pluralities of HMIX, FRPE, and TABR2 had no plots available for analyses and are not expected to have sufficient 
MPB host to warrant a hazard above 0 so are not included in this table. 
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| APPENDIX B (cont) 































































































































































MOUNTAIN PINE BEETLE 
| In WESTERN WHITE PINE** 
HAZARD RATING** SAMPLE SIZE** 
| Dominance Size Class (inches DBH | Dominance Size Class (inches DBH) 
' Group 60% | 0.0-4.9 | 05.0-09.9 | 10.0-14.9 | 15.0+ Group 60% | 0.0-4.9 | 05.0-09.9 | 10.0-14.9 | 15.0+ 
Plurality! Total Total Total Total Plurality! Total Total Total Total 
ABGR 0 1 0 ABGR 10 34 65 
ABLA 0 0 0 0 ABLA 83 88 
BEPA | BEPA 
| CELE3 CELE3 
| IMIX l 0 0 0 IMIX 
JUNIP JUNIP 
LALY LALY 
LAOC 0 0 LAOC 
PIAL 0 PIAL 
| PICO 0 0 0 PICO 
| PIEN 0 0 PIEN 
PIFL2 PIFL2 
PIMO3 PIMO3 
PIPO 0 0 PIPO 
| POPUL POPUL 
| POTRS POTRS 
PSME 0 0 0 0 PSME 
| THPL 0 THPL 
| TMIX 1 1 0 0 TMIX 
TSHE 1 TSHE 
TSME 0 TSME 13 














*only applicable west of Continental Divide where white pine occurs 


** Color coding was used to indicate sample sizes; estimates based on >10 plots has no color; on 1-9 plots are gold, and on 0 plots are orange. 
| Hazard ratings based on <10 plots were determined using expert opinion. For calls based on 1-9 plots, changes to predicted hazard based on 
professional opinion are indicated as the suggested hazard with predicted in () following [e.g. “1 (not 0)”]. 


‘Dominance Group 60% Pluralities of HMIX, FRPE, and TABR2 had no plots available for analyses and are not expected to have sufficient MPB host to 
warrant a hazard above 0 so are not included in this table. 
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APPENDIX C: Alterations to model for MPB in whitekark and limber pines based on | 
Continental Divide 


MOUNTAIN PINE BEETLE | 
In WHITE BARK and LIMBER PINES | | 


HAZARD RATING*? 


Dominance Size Class (inches DBH 


Group 60% 
0.0-4.9 Total eee 0-09.9 Total 10.0-14.9 Total ens 0+ Total 


Plurali 
ABGR 
== 1 (E22 Sar any eS 
1 (W=0) 



























BEPA 
CELES 


PICO 


PSME 
THPL 


TMIX 
TSHE 
TSME 


*Color coding was used to indicate sample sizes; estimates based on >10 plots has no color; on 1-9 plots are gold, and on 0 plots are orange. 
Hazard ratings based on <10 plots were determined using expert opinion. For calls based on 1-9 plots, changes to predicted hazard based on 
professional opinion are indicated as the suggested hazard with predicted in () following [e.g. “1 (not 0)”]. - | 




















































'Dominance Group 60% Pluralities of HMIX, FRPE, and had no plots available for analyses and are not expected to have sufficient MPB | 
host to warrant a hazard above 0 so are not included in this table. | 


Hazard rating changes (GREEN FONT) from Regional model (BLACK FONT) for East (E) and (W) of the Continental Divide; values in 
green that are “suggested” have fewer than 10 plots for decision making. 
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APPENDIX C (cont.) 
MOUNTAIN PINE BEETLE 
In WHITE BARK and LIMBER PINES 


Dominance Dominance 
: Total Total Total Total Plurality! Total Total Total Total 
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*Color coding was used to indicate sample sizes; estimates based on 

>10 plots has no color; on 1-9 plots are gold, and on 0 plots are orange. Hazard ratings based on <10 plots were determined using expert 
opinion. For calls based on 1-9 plots, changes to predicted hazard based on professional opinion are indicated as the suggested hazard with 
predicted in () following [e.g. “1 (not 0)”]. 


'Dominance Group 60% Pluralities of HMIX, FRPE, and TABR2 had no plots available for analyses and are not expected to have sufficient 
MPB host to warrant a hazard above 0 so are not included in this table. 
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APPENDIX D: Exploring forest level adjustments 


The Regional level MPB hazard ratings for ALL PINE, LPP, PP, WBkLmP and WWP may be re- 
evaluated at the Forest Level if the REGION model has >10% error in a hazard level of interest (Table 9). 
The following is an example of some of the analyses and changes that were explored to create what we 
hope are more accurate depictions of MPB stand susceptibility. In these examples the ALL PINE model 
was used although any one of the host-specific models (LPP, PP, WBkP-LmP, and WWP) may be 
similarly evaluated and altered. This process should be done cooperatively with Forest personnel 
involved in landscape analyses and an FHP entomologist familiar with this modeling. 


Model adjustments based on Forest-specific R1-SDB plots — the Gallatin NF example 


Table A shows the results of analysis using only FIA plots found on the Gallatin NF. Of the nine possible 
changes suggested by the Gallatin-only data, most are in Dom60 vegetation types that have been 
identified as having a bimodal distribution Regionally or at least a wide range of hazard levels but only 
three are based on 10 or more plots. 


Changes to the Regional model hazard rating should be considered carefully and justifications recorded. 
In this Gallatin NF case, we decided up front that three suggested changes were based on too few of plots, 
did not come from Dom60 types likely to have highly variable data, and were not helpful in highlighting 
areas most likely to have high susceptibility to MPB (i.e. changes in PIAL, PICO and PIFL2). Table B 
provides the justifications used for the other five changes in ALL PINE hazard rating specific to the 
Gallatin NF. 


Making the five suggested changes results in a better overall model of hazard rating for the Forest. Table 
C shows the Gallatin ALL PINE model improvement over the Regional ALL PINE model (also shown in 
Table 9). A closer look at plot-level mismatch in hazard rating indicates that the ratings given to FIA 
ground plots using Randal and others (201 1b) were correctly predicted by our model using R1-ExVeg 
characteristics about 60% of the time. However, this does not include the expected decrease in overall 
error if hazard levels VERY LOW (0) and LOW (1) are combined. A few instances where 1) mismatches 
were overestimated more than one level in the model (-2) or 2) mismatches were consistently 
underestimated by the model >25% of the time are highlighted. Further exploration of these mismatch 
errors may be warranted depending on the modelling goals or the composition of the analyses area within 
the Gallatin NF. 
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TABLE A: GALLATIN NF HAZARD RATING 
for ALL PINE HOST SPECIES 










































HAZARD RATING* : SAMPLE SIZE* 
Dominance Size Class (inches DBH) Dominance Size Class (inches DBH) —__ | 
Group 60% [| 0.0-4.9 | 05.0-09.9 | 10.0-14.9 | 15.0+ Group 60% | 0.0-4.9 | 05.0-09.9 | 10.0-14.9 | 15.0+ | 
Plurality! Total Total Total Total Plurality! Total Total Total Total | 
R eet eo ig ie a oe. a ea 
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POPS eo GS oe 

| CELE3 Pee mee oa iin 
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*Purple numbers indicate Regional ALL PINE hazard ratings. Larger 

red numbers indicated changes suggested by Gallatin NF-only plot analyses that were accepted and changed; blue numbers indicate changes 
suggested but not changed. Color coding of cells was used to indicate sample sizes; estimates based on >10 plots has no color; on 1-9 plots 
are gold, and on 0 plots are orange. If the Dom60 group did not appear in the Gallatin FIA plot dataset it is colored gray. Hazard ratings 
based on <10 plots were determined using expert opinion. For calls based on 1-9 plots, changes to predicted hazard based on professional 


opinion are indicated as the suggested hazard with predicted in () following [e.g. “1 (not 0)”]. 


'Dominance Group 60% Pluralities of HMIX, FRPE, and TABR2 had no plots available for analyses and are not expected to have sufficient 
MPB host to warrant a hazard above 0 so are not included in this table. 


TABLE B.: justification for the five changes recommended for the Gallatin NF, ALL PINE hazard rating in R1- 
VMAP 























Reasoning for recommended change 


a move from 0 to 3 is too large based on one plot; ABLA on the Gallatin 
NF is not typically considered high hazard to MPB; changing to a 1 to be 


Dom60 
vegetation | Size Change 
ae we made 
more in line with the rest of ABLA type is reasonable although O and 1 


change from 
Otol 
hazard levels may be combined so may be unnecessary 
IMIX 05.0- change from | IMIX is a vegetation type that by definition is going to have wide 
ee Yeas 2to3 variability in composition over the Region; IMIX on Gallatin appears to 
have more pine host than IMIX averaged over the Region 
10.0- change from | IMIX is a vegetation type that by definition is going to have wide 
14.9 2 to3 variability in composition over the Region; IMIX on Gallatin appears to 
have more pine host than IMIX averaged over the Region 
15.0+ change from | IMIX is a vegetation type that by definition is going to have wide 
2 to 3 variability in composition over the Region; IMIX on Gallatin appears to 
have more pine host than IMIX averaged over the Region 
PSME PSME is strongly bimodal vegetation type; for the Gallatin a 0 hazard is 
eA. | 14.9 2 to 0 more appropriate than a 2 hazard 
TMIX change from | TMIX is a strongly bimodal vegetation type that, over the Region, had an 
average hazard rating of 1; for the Gallatin a 2 hazard rating is more 
appropriate 


TABLE C: Comparison of Y%carea predicted to be in each of three hazard rating classes based on plot-level hazard 
estimates and R1-ExVeg classification hazard estimates using the ALL PINE hazard model and FIA periodic plot 
data (1993-2002) for the Northern Region (REGIONAL) or Gallatin National Forest (FOREST). | 
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TABLE D: A closer evaluation of the mismatch between the hazard ratings given to FIA plots located on the 

Gallatin NF (“reality”) and the hazard een to me R1-VMAP category. of that plot (“modelled”). 
modelled as rer | 

higher hazard 
than 9 a) a 
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31% | 38% | 














| PSME still contains a wide range | 
| of hazard levels and may need 


|__| 100% 
| further evaluation depending on 


| sae 37) 46 | modelling goals 

| | Review of plot data shows TMIX 
| | | hazard plots were mis- 

| TMIX 5% | 63%| 21%| | 38 | 


10-14.9” is where many high 
| identified; this may warrant 
[ Grand 
| Total | 5%] 22% 
























further evaluation depending on | 
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Model adjustments based on R1-VMap polygons - the Lolo NF project area example 


Reassessing hazard ratings for areas smaller than a National Forest should not be done by reanalyzing FIA plots for 
that area; there will likely be too few to base changes on. In some cases both the Regional models and the Forest- 
level models will have vegetation types (Dom60) that contain more error than is acceptable for the smaller analyses 
area. One possible way to further refine hazard for these Dom60 types of concern is to examine the R1-VMap data 
for that analyses area. More specifically, each polygon contains information not only on the Dom60 vegetation 
class but the more detailed DG6040 which can be explored. 


An analysis similar to what was done on the Gallatin was completed for the Lolo National Forest for the ALL 
PINE, LPP, and PP hazard models. Most potential changes in hazard rating at the Forest level appeared 
inconsequential (e.g. between ratings of VERY LOW [0] and LOW [1]), had few plots for basis (<10), and did not 
reflect changes potentially justified by bark beetle stand preferences. The few potential alterations we thought 
might be worth exploring included: 


ALL PINE MODEL 
-' ~PSME)15*+ possible increase to hazard of 2 (from 1) 
- TMIX 10-14.9” and 15”+ issues of bimodal data (ratings of 1 versus 2) 
- IMIX sizes >4.9” no hazard rating changes indicated but most plots on the Lolo with hazards 
of 3 that were rated lower by the model were in this group 
PP MODEL 
=) PIPOW0-l 55 possible change to hazard of 3 (from 2) but only 3 plots 
- IMIX sizes >9.9” no hazard rating changes indicated but strong bimodal distribution of data 
between hazard ratings 0 and 2 
LPP MODEL 


- No changes indicated but noticed high bimodal tendency remained in IMIX 10-14.9” 


Before making changes to the hazard ratings applied to R1- VMap we wanted to know if any of these Dom60 
vegetation types were common or important to the analysis area. We discovered that IMIX was both common and 
very likely to have pines as part of the species composition (versus other shade intolerant tree species that are not 
typical MPB hosts). 


Within the analysis area our GIS coverage indicated 18,373 acres of IMIX type (Dom6). Because IMIX is a 
heterogeneous category where no single species makes up 60% or more of the basal area, greater detail on presence 
of pine can be gotten from DG6040 classifications. Table E shows this additional breakdown within the IMIX 
(Dom60) type: 













Table E. Distribution of IMIX Dom60 vegetation type among Dom6040 vegetation types 


Dom6040 type %of Project Area 


IMIX No one tree species has 40% dominance or greater 


LAOC has 40-59% dominance; the remainder is mixed, jase | see 
LAOC-IMIX shade intolerant species, none of which constitute 40% or 2,954 16.1% 

more of the basal area 
PICO-IMIX PICO has 40-59% dominance (remainder as above) 
PIPO-IMIX PIPO has 40-59% dominance (remainder as above) 
PSME-IMIX PSME has 40-59% dominance (remainder as above) 


Grand Total (Dom60IMIXtype) | 18373 | 


: 
: 


Description 
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If within the IMIX type, only the DG6040 types of PICO-IMIX and PIPO-IMIX types have the potential for MPB 
activity, we determined that approximately 30% of the IMIX type within the project area could have a HIGH 
hazard rating (15.7%+14.4%). We then broke down the PICO-IMIX and PIPO-IMIX types further by their canopy 
cover and assumed that canopy covers less than 40% would constitute stands of densities too low to be highly 
susceptible to MPB (Table F). This summary suggests that only 14% of the IMIX (Dom60) type has 40-59% pine 
(PICO or PIPO basal area), and over 40% canopy cover. If we assumed only the larger diameters (>9.9’’) and 
greatest density stands (>60% canopy) fully meet high hazard thresholds, the area of concern for MPB activity 
drops to around 3.4% of IMIX in the project area. 





Table F. The breakdown of Dom60 vegetation type IMIX into subcategories most likely to have high 
susceptibility to mountain pine beetle. This includes the two pine-dominated DG6040 vegetation types 
PICO-IMIX and PIPO-IMIX, the denser stands with CANOPY COVER of 40% or greater, and the SIZE 
CLASS of 5”DBH or greater. Red type denotes what is likely the highest hazard areas. However, 
CANOPY COVER and SIZE CLASS are descriptive of the overall stand and not of the pine host type 
separately. 























% of DG6040 %ofDG6040 %of % of 

types PICO- types PICO- Dom60 Dom60 
CANOPY TREE ACRES in IMIX + PIPO IMIX + PIPO type type 

Dom60 DG6040 COVER SIZE PROJECT IMIX IMIX IMIX IMIX 












IMIX | PICO-IMIX | 40 to 60% 71.5% 10.0% 
560% | 
PIPO-IMIX | 40 to 60% 28.5% 4.0% 
om 
ke. ae eee 
Grand Total 2573 100% 14% 





Evaluation of Rl-VMap composition in the project area suggests that the Dom60 IMIX vegetation type doesn’t 
contain many stands considered HIGH (3) hazard to MPB activity. However, if those analyzing the project area 
want to be liberal about identifying stands susceptible to MPB, they may choose to highlight any R1-Veg types that 
have a significant possibility of having HIGH (3) hazard. To make this change to R1-VMap they would adjust 
MODERATE (2) hazard rating to HIGH (3) in the attribute table for the IMIX (Dom6) type for all but the 0-4.9” 
size class. Table G is an example of the justification for these changes; information that should be reported with 
any changes from the Regional model. 
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Table G. Various traits of the IMIX (Dom60) vegetation type for the Northern Region, for the Lolo National 
Forest, and for a given project area. This information can be used to determine if specific changes should be 
made to ALL PINE hazard ratings within the project areas. 
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Around % of the FIA plots on the Lolo NF that fall into these 
05.0- 6058 21% 20r3 three IMIX-size classes were given HIGH hazard (3) ratings. 
09.9 Thus a liberal effort to identify areas in VMAP with any 
likelihood of high hazard could change ratings from MODERATE 
10.0- 6846 5 5 31 26% aOna (2) to HIGH (3). We also noted that these 3 categories have 
14.9 
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~26% of the treed acres in the project area so they are 
important. However, when we look at all IMIX within the 
Project area in each size class, and choose only those stands 
with pine (PIPO or PICO) constituting at least 40% of plot, and 
canopy cover >=40%, we find only 6%, 6%, and 2% of all IMIX 
acres likely warrant a HIGH rating, respectively (by size class). 
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